Towards optimal softening in three-dimensionalN-body codes - I. Minimizing the force error
نویسندگان
چکیده
منابع مشابه
Towards optimal softening in 3 D N - body codes : I . Minimizing the force error
In N -body simulations of collisionless stellar systems, the forces are softened to reduce the large fluctuations due to shot noise. Softening essentially modifies the law of gravity at r = |xi − xj | smaller than some softening length ǫ. Therefore, the softened forces are increasingly biased for ever larger ǫ, and there is some optimal ǫ which yields the best compromise between reducing the fl...
متن کاملAn energy-conserving formalism for adaptive gravitational force softening in SPH and N−body codes
In this paper we describe an adaptive softening length formalism for collisionless N−body and self-gravitating Smoothed Particle Hydrodynamics (SPH) calculations which conserves momentum and energy exactly. This means that spatially variable softening lengths can be used in N−body calculations without secular increases in energy. The formalism requires the calculation of a small additional term...
متن کاملOptimal Softening for N -body Simulations
We propose a strategy for determining the optimal softening length in N -body simulations. This strategy ensures that the ensemble-average acceleration error is minimal at a small radius r0, and that the error never exceeds the optimal value beyond r0. Analytic formulae are derived for calculating the acceleration error for any halo profile or softened gravity so that the optimal softening leng...
متن کاملOptimal Softening for N -body Halo Simulations
We propose to determine the optimal softening length in N -body halo simulations by minimizing the ensemble-average acceleration error at a small radius r0. This strategy ensures that the error never exceeds the optimal value beyond r0. Furthermore, we derive semi-analytic formulae for calculating the acceleration error due to the discreteness of particles and softened gravity, which are valida...
متن کاملTowards Optimal Parallel PM N-body Codes: PMFAST
We present a new parallel PM N-body code named PMFAST that is freely available to the public. PMFAST is based on a two-level mesh gravity solver where the gravitational forces are separated into long and short range components. The decomposition scheme minimizes communication costs and allows tolerance for slow networks. The code approaches optimality in several dimensions. The force computatio...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
ژورنال
عنوان ژورنال: Monthly Notices of the Royal Astronomical Society
سال: 2001
ISSN: 0035-8711,1365-2966
DOI: 10.1046/j.1365-8711.2001.04237.x